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The Formation of Distributed and Clustered Stars 
in Molecular Clouds 



1.1 Introduction 

During the last two decades, the focus of star formation research has shifted 
from understanding the collapse of a single dense core into a star to studying 
the formation hundreds to thousands of stars in molecular clouds. In this 
chapter, we overview recent observational and theoretical progress toward 
understanding star formation on the scale of molecular clouds arid coin - 
plexes, i.e the macrophysics of star formation ( McKee and Ostrikeil . 20071 ). 



We begin with an overview of recent surveys of young stellar objects (YSOs) 
in molecular clouds and embedded clusters, and we outline an emerging pic- 
ture of cluster formation. We then discuss the role of turbulence to both 
support clouds and create dense, gravitationally unstable structures, with 
an emphasis on the role of magnetic fields (in the case of distributed stars) 
and feedback (in the case of clusters) to slow turbulent decay and medi- 
ate the rate and density of star formation. The discussion is followed by 
an overview of how gravity and turbulence may produce observed scaling 
laws for the properties of molecular clouds, stars and star clusters, and how 
the observed, low star formation rate may result from self regulated star 
formation. We end with some concluding remarks, including a number of 
questions to be addressed by future observations and simulations. 



1.2 Observations of Clustered and Distributed Populations in 

Molecular Clouds 

Our knowledge of the distribution and kinematics of young stars, protostars 
and dense cores in molecular clouds is being rapidly improved by wide field 
observations at X - ray, optical, infrared , and (sub)millimeter wavelengths 



(j Allen et al.l . 120071 : iFeigelson et al.l . 120071 ) . In this section we present multi 
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wavelength portraits of some of the most active star forming regions near 
our Sun and discuss some imphcations of these results. 



1.2.1 Molecular Cloud Surveys 

The Spitzer space telescope, with its unparalleled ability to detect IR-excesses 
from YSOs with disk and envelopes over wide fields, is now producing the 
most complete censuses to date of young stellar objects (YSOs) in nearby 
molecular clouds. One example is the survey of the Orion A cloud, the 
most active star forming cloud within 450 pc of the Sun (Megeath et al., in 
prep). The observed distribution of YSOs in this cloud exhibits structure 
on a range of spatial scales, tracing both the overall filamentary morphology 
of the clouds and smaller filamentary and clumpy sub-structures within the 
clouds, as shown in Fig. 11.11 This figure also shows a range of observed 
stellar densities; YSOs can be found in large clusters containing OB stars, 
in smaller groups and in relative isolation. 

Is there a preferred environment for star formation in giant molecular 
clouds: do most low mass stars form in rare large clusters, the more numer- 
ous small groups, or in a distributed population of more isolated stars? To 
some extent, this depends on the criteria used to distinguish isolated and 
clustered stars; as shown in Orion A, the clusters are often not distinct ob- 
jects but belong to an extended distribution of stars (Fig. ll.ip . One approach 
is to adopt a critical surface density for clustered stars, although there is no 
clear motivation for a particular threshold density as demonstrated by the 
continuous distribution of YSO surface densities shown in Fig. II. li In the 
Orion cloud complex (Orion A and B), Megeath et al. (in prep) decomposed 
the observed distribution of YSOs with IR-excesses into groups of 10 or more 
sources where each member is in a contiguous region with a local surface 
density of 10 stars per pc~^ or higher. They find that 44% are in the 1000 
member Orion Nebula Cluster (ONC), 18% are in seven clusters with 30 to 
100 YSOs and 9% are in twelve groups of 30 to 10 YSOs. They find that 
28% of the YSOs with IR-excesses (639 objects in t otal) form a distr ibuted 
population outside the groups and clusters (also see lAllen et al.l . l2007l l. The 
number of YSOs in the ONC is underestimated due to incompleteness, and 
consequently the actual fraction of stars in the distributed population is 
somewhat lower (Megeath et al. in prep). Despite this bias, the surveys are 
showing that even in giant molecular clouds (GMCs) forming massive stars, 
a substantial number of YSOs are found in relative isolation; the observed 
range of densities and environments must be explained by models of cloud 
fragmentation and evolution. 
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Fig. 1.1. Left: The distribution of nearest neighbor densities {10/{'kE^q), where 
i?io is the distance to the 10th nearest neighbor) for the Orion A GMC and the 
Taurus dark cloud complex (K. Luhman, priv. comm.). The dotted line gives the 
adopted threshold density for separating distributed stars from stars in groups and 
clusters. Right: The distribution of Spitzer identified young stellar objects with 
infrared excesses in the Orion molecular clouds. The gray line gives the outline of 
the survey fields, including the large hockey stick shaped field and three reference 
fields. The black lines is the Ay = 3 contour determined from an extinction map 
constructed from the 2MASS point source catalog ( R. Gutermuth , priv . comm.). 
The star symbols give the positions of OB stars from iBrown et all (|1994( ). 



1.2.2 A Gallery of Embedded Clusters 

In Figs. 11.21 11.31 11.41 & 11.51 we show the latest compilations of YSOs in 
four of the nearest embedded (or partially embedded) clusters to the Sun. 
The most striking result is the diversity of configurations apparent in the 
clusters. The dense gas in the Serpens cluster is concentrated in a 0.5 pc long 
filament divided into two main clumps; the protostars are primarily found 
in the two clumps. In contrast, the pre-main sequence stars are more evenly 
distributed over a 1 pc diameter region with only a modest concentration 
in the central region; this may be the result of dynamical evolution due 
to gravitational interactions between YSOs and due to the non-spherically 
symmetric cloud potential. The dense gas in NGC 1333 forms a network 
of clumps and filaments spread over a region more than 1 pc in diameter. 
Both the protostars and pre-main sequence stars follow this structure, with 
the protostars concentrated in the dense gas and pre-main sequence stars 
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Protostars Pre — Main Sequence Stars 
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Fig. 1.2. The distributio n of young stel lar objects in the Serpens cloud core cluster 
at a distance of 260 pc (j Winston et al .. 2007). The contours show the dense gas 
detected in an 850 fim SCUBA map feavis et al.l . I2000f ). The left panel shows the 
distribution of protostars separated into flat spectrum (gray markers) and Class I/O 
(black markers) sources. The right panel show the pre-main sequence stars with 
disks (black markers) and without disks (gray markers). The pr e-main sequence 
stars with out disks were identifi ed by their elevated X-ray emission (jCiardino et al.L 
[2003) ; see I Winston et al.l (|2007t) for a discussion of the completeness the displayed 
source selection. Note that the SCUBA map and Chandra images do not cover the 
entire displayed region. 



found in the immediate vicinity of the dense gas; this implies little dynamical 
evolution in this region. Unlike the previous two clusters, the IC348 cluster 
shows a circularly symmetric and centrally condensed distribution. The 
molecular gas has been largely cleared from this region; the protostars are 
concentrated in a filamentary gas structure on the edge of the cluster. 

The massive Orion Nebula Cluster (ONC) shows a particularly complex 
morphology (Figs. 11.11 fc Fig. II. 5p . In the center of the ONC is the densest 
known region of young stars in the nearest 500 pc. This central conden- 
sation is elongated and aligned with the axis of the filamentary molecular 
cloud. To the north of the central condensation a filamentary distribution of 
protostars follows the molecular cloud. It is surrounded by a more extended 
distribution of pre-main sequence stars. To the south, the molecular cloud 
appears to have been partially swept up into a shell by the massive stars 
in the center of the ONC. A more extended distribution of YSOs in this 
regions fills the shell, indicating that these stars probably originated in the 
expanding shell. 



1.2 Observations of Clustered and Distributed Populations in Molecular Clouds^ 



Protostors Pre — Main Sequence Stars 




3:30:00.00 3:29:00.00 3:28:00.00 3:30:00.00 3:29:00.00 3:28:00.00 

RA (2000) RA (2000) 



Fig. 1.3. The distribution of YSOs in the NGC 1333 cluster at a distance of 
250 pc. The contours s how the dense gas traced in the 850 SCUBA map 
(jSandell and Kne3 . 1200 if ). The left panel shows the distribution of protostars 
(Gutermuth et al., 2007); the right panel shows the distribution of pre-main se- 
quence stars with disks (black markers; Gutermuth et al., 2007); a nd diskless pre- 
main sequence stars with elevated X-ray emission (gray markers; iGetman et al.l 
|2002[ ). Note that the Chandra survey of Getman et al. does not cover the entire 
displayed field. 



1.2.3 The Distribution of Protostars in Embedded Clusters 

The distribution and spacing of protostars is an important constraint on 
the physics of fragmentation and the potential for subsequent interactions 
between protostars. As shown previously, the protostars trace the clumpy 
and filamentary distribution of the dense molecular gas. The protostars in 
embedded clusters are often closely packed, with med ian projected nearest 
neighbor spacings ra nging from 5000 AU to 20000 AU (jWinston et al.l . 120071 : 



Muench et al.l . 120071 ). In a few cases, peaks are fo und in the distribution of 



projected nearest neighbor distan ces for protostars (iGutermuth. R. A. et al 



20071 : iTeixeira. P. S. et al.l . I2OO6I ) and for dense cores (jStanke et al.1 . boOQ ) 



These preferred separations for the protostars are within a factor of three 
of the local Jeans length, suggesting that the protostars may result from 
gravitational fragmentation (Teixeira et al., 2006). The detection of ex- 
tremely dense groups of protostars with ~ 7 objects in regions ~ 10000 AU 
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Fig. 1.4. The distribution of YSOs in the IC 348 Cluster. The contours show the 
dense gas as traced by the 850 /im SCUBA map (J. DiFrancesco, priv. comm.). 
The left panels show the distribution of protostars (black markers); the right panel 
shows the distribution of pre-main sequence stars \ vith (black markers) a nd without 
(gray markers) disks. The sample was taken from iMuench et aTl (|2007t ). 



diameter sugge s t that hierarchical fragme ntation operates in some cases 



(jWinston et al.1 . l2007l : iTeixeira et alJ . l2007l ) 



Can the protostars in the observed dense groups interact? In the Ser- 
pens region the close spacing of protostars implies that the volumes from 
which these protostar s accrete are often densely packed, if not overlapping 
, 20071 ) ■ In a study of a dense group of proto stars in Serpens , 
( 20071 ) argued that the high gas density (jOlmi and Testil . 

as measured in 



(Winston et al. 



Winston et al 

2002h . and the subvirial RMS velocity of the protostars 



the B IMA N2H"'" detections of the protostellar cores; IWilliams and Mvera . 
2OO0I ) are consistent with the protostars accreting a significant portion of 
their mass through Bondi-Hoyle accretion. Bondi-Hoyle accretion in such a 
closely packed group may lead to the cor npetitive accretion predic ted in some 



numerical models of turbulent clouds ( Bonnell and Bate . 20061 ). Further- 



more, if such groups st ay bound for 400,000 years (the estimated duration 



of the protostellar phase iHatchell et all 120071 ). each large (1000 AU) proto- 



stellar envelope would typically experience one collision. Further studies are 
needed to find direct evidence for interactions and to ascertain the fraction 
of stars that form in dense groups where the potential for interactions exists. 
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Fig. 1.5. The distribution of YSOs in the Orion Nebula Cluster. The contours show 
an extinction map constructed from the 2MASS PSC (R. Gutermuth, priv. comm.). 
Left panel: The large black circles are protostars identified in the Spitzer data 
(Megeath et al., in prep.); no protostars are identified in the center of the Orion 
nebula where the 24 /im data is saturated. The small blac k dots are the detecte d 
X-ray sources from the Chandra Orion Ultradeep Project (|Feigelson "eralll200l : 
this survey covered a diamond shaped region defined by the extent of the X-ray 
sources. Right Panel the distribution of Spitzer identified pre-main sequence stars 
with disks (black markers; Megeath et al. in prep) and pre-main sequence stars 
without disks (gray markers) identified by variability ([Carpenter et al.l . l200lD . 



1.2.4 The Evolution of Embedded Clusters: An Observational 

Perspective 

A common trait of embedded clusters is a lack of circul ar symmetry: they 



are o f ten elongated and in many cases show substructure (jGutermuth et al 



2OO5I : lAllen et al.l . 120071 ) . This indicates that clusters retain the imprint of 



the filamentary, clumpy structure of the molecular clouds from which they 
formed; consequently they do not appear to be relaxed, virialized clusters 
of stars. This is supported by observations of the kinematics of protostel- 



observation of the NGC1333 and NGC2264 clusters; 



Peretto et al.1 . l2007l ) 



clusters; 


Walsh et al. 


1 


2007; 


Furesz et al. 


( 


2007 


) find that 
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the velocities of the pre-main sequence stars largely follow the overall veloc- 
ity gradient in the gas. This shows that the more evolved pre-main sequence 
stars also retain the primordial velocity structure of their parental cloud. 

Observations of 3-5 Myr clusters sh ow that the paren tal gas is largely 
dispersed (see discussion of examples in lAUen et al.l . 120071 ): at this time the 
cluster may become unbound since the majority of the binding ma ss has 
been ejected from the cluster ( Lada et al. . 1984 : Kroupa et al. . 200ll ). The 
question is whether this timescale allows for the cluster to relax. Adopting 
radii of 0.5 pc for the clusters (encompassing the dense central regions of 
the clusters) and using the estimated number of stars (420, 156, 182 and 
1000 stars for IC 348, NGC 1333, Serpens and ONC), we estimate relax- 
ation times of 2.4, 1.8, 1.8 and 6 Myr, respectively. The actual relaxation 
times will be longer since the mass in the clusters is dominated by the gas 



(jAdams and Mvera . l200ll ). This suggests that the relaxation time is ap- 



proximately equal to or less than the gas dispersal time. Star formation 
continues in the residual cloud as the gas is being di s perse d, as shown in 
IC 348 and the ONC (Fig. O & [LSD . iMuench et all (|2007l l estimate that 
the current star formation rate in IC 348 is consistent with a constant star 
formation rate over the entire lifetime of the cluster. Thus, the timescale for 



star formation, gas dispersal and relaxation are similar (cf. also lAllen et al 



20071 ). 



A picture is emerging where clusters form in filamentary, clumpy clouds. 
During the first few million years of their evolution, clusters form stars over 
regions 1 pc in diameter. During this time, the surrounding gas is being 
dispersed by winds and radiation, eventually leading to the cessation of star 
formation in 3-5 Myr. Although the stars may form with a low velocity 
dispersion, at least relative to local structures in the ga s, collapse and large 
scale velocity gradi ents in clouds may stir up the cluster (jPeretto et al.l . 120071 : 
Furesz et alll2007l ). In extr eme cases, this may lea d to violent relaxation, as 
suggested for the ONC by Feigelson et al. ( 20051 ). Dynamical interactions 
between stars and between the stars and the cloud potential will scatter 
stars; however, given that the cloud mass (which typically embodies 3/4 of 
the total mass) is being ejected on a timescale similar to the relaxation time, 
it is unlikely that most embedded clusters will be able to achieve a relaxed 
equilibrium. One counter example may be IC 348, where the circularly 
symmetric, c entrally condensed c onfiguration of this cluster suggests that 
it is relaxed (jMuench et al.l . 120071): however, circular syin metry could also 
result from dynamical expansion (IGutermuth et al.1 . [2OO5I 1 . N-body simula- 
tions are needed to fully understand what the observed cluster morphologies 
imply about the dynamic state of the cluster. 
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1.3 Local Theory of Distributed and Clustered Star Formation 

Stars form in turbulent, magnetized, clouds. The relative importance of 
magnetic fields and turbulence in controlling star formation is a matter 
of debate. Early quantitative models have concentrated on the forma- 
tion and evolution of individual (quiescent, low-mass) cores and the role 



of magnetic fields ( Shu et al. . 1987 : Mouschovias and Ciolek . 19991 ). More 



recent studies have concentrated on the role of turbulence in cloud dynarn - 
ics and core formation, as r eviewed in, e.g.. iMac Low and KlessenI (120041): 



Elmegreen and Scald (j2004l ): lBallesteros-Paredes et al.l (j2007l ). and lMcKee and Qstriker 



(j2007l ). Ultimately, the debate can only be settled by direct measurement of 
the flux-to-ma ss ratio for the bulk of the molecular gas, which is currently 
not available ( Crutcher . 19991 ). In the absence of such measurements, we 
h ave to rely on indirect evidence and theoretical arguments. 

McKeel (| 19991 ) argued that self-gravitating GMCs should be magnetically 
supercritical by a factor of ~ 2 as a whole, if turbulen ce and orde red ma; 



netic fields provide comparable cloud support (e.g., iNovak et al.l . |200 



Even in such a globally (moderately) supercritical GMC, the magnetic field 
may still play a key role in regulating star formation, because the formation 
of a typical star or small stellar groups involves only a small sub-piece of 
the cloud that contains a few to tens of solar masses. Such a region can be 
locally subcritical, with a mass less than the magnetic critical mass, w hich 
can be hundreds of solar masses for typical parameters ( McKee . 19991 ). un- 
less (1) a large fraction (perhaps ~ 1/2 or more) of the material along a flux 
tube has collected in the region, or (2) the region is part of a substructure 
(i.e., a dense clump) that is already highly supercritical to begin with. In 
the former case, it would take a long time (perhaps more than 10 Myrs) 
for the required mass accumulation to occur in turbulent GMCs of tens of 
parsecs in size and a few km/s in rms speed, if it happens at all. If not, 
locally supercritical sub-pieces of stellar masses must be created by some 
other means, most likely through ambipolar diffusion. In the latter case, 
ambipolar diffusion is not expected to play a decisive role in star formation. 
The two cases may corr espond to, respe ctively, distributed and clustered 



modes of star formation ( Shu et al. . 198?! ) 



1.3.1 Ambipolar Diffusion and Distributed Star Formation 

The best studied region of distributed star formation is the Taurus molec- 
ular clouds. Although it is not clear whether these clouds are magnetically 
supercritical or subcritical as a whole, the dynamics of the more diffuse re- 
gions is probably magnetically dominated. The best evidence comes from 
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thin strands of ^^CO emission that are ahgned with the directions of the lo- 
cal magnetic field (Heyer et al. 2008, preprint), where the field is apparently 
strong enough to induce a measurable difference between the turbulent ve- 
locities along and perpendicular to the field direction. The filamentary mor- 
phology is strikingly similar t o that observed in the nearby Riegel-Crutcher 
HI cloud, mapped recently by McClure-Griffiths et al. ( 20061 ) using 21cm ab- 
sorption against the strong continuum emission towards the galactic center. 
Its filaments are also along the directi ons of the local magnetic fi eld, which 
has an estimated strength of ~ 30 AiG. iKazes and Crutcheil dlQSd ) measured 
the line-of-sight field strength at a nearby location, and found Bios ~ 18 /iG, 
which is consistent with the above estimate if a correction of ~ 2 is applied 
for projection effect. The inferred strong magnetization of this cloud may 
not be too surprising, in view of the result that cold neutral HI structure s 
are strongly magnetically subcritical in general ( Heiles and Troland . 20051 ). 
If diffuse molecular clouds such as the Taurus clouds are formed out of 
such HI gas, it is not difficult to imagine that at least some of them will 
be subcritical as well. If the bulk of a molecular cloud is indeed magnet- 
ically subcritical, then the well-known low efficiency of star formation can 
be naturally explained: the formation of dense, star-forming, magnetically 
supercritical cores is regulated by ambipolar diffusion, which is generally a 
slow process, with a tir ae scale an order of magnitude (or more) longer t han 
the local free fall time ( Shu et al. . 1987 : Mouschovias and Ciolek . 19991 ). 

For relatively diffuse molecular gas with an Ay of a few, the problem may 
be the opposite: the ambipolar diffusion time scale is too long to allow for 
significant st ar formation i n a reasonable time due to io nization from UV 
background (McKee, 19891 : Mvers and Khersonskv . 19951 ). In order to form 
stars in a reasonable time scale, the rate of ambip olar diffusion must be 
enhanced (jZweibell . |2002| : iFatuzzo and Adama . |2002| ) . This is where turbu- 
lence can help greatly. Supersonic turbulence naturally creates dense regions 
where the background UV photons are shielded, and where the gradient in 
magnetic field is large; thereby accelerating ambipolar diffusion. The strong 
magnetic field, on the other hand, prevents too large a fraction of the cloud 
mass being converted into stars in a turbulence crossing time. In this hybrid 
scenario of distributed star formation in relatively diffuse molecular clouds, 
both magnetic fields and turbulence play crucial roles: the star formation 
is accelerated by turbulent compression but r e gulat ed by magnetic fields 
through ambipolar diffusion (|Li and Nakamural. 120041'). 

The above scenario is illustrated in Nakamura and Li ( 20051 ). using 2D 
simulations of a sheet-like, magnetically subcritical cloud (with an initial 
dimensionless flux-to-mass ratio F = 1.2 everywhere). It is stirred with a 
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(compressive) supersonic turbulence of rms Mach number A4 = W at t = 0. 
The strong initial magnetic field prevents the bulk of the strongly shocked 
material from collapsing promptly. The shocked material is altered perma- 
nently, however; its flux-to-mass ratio T is reduced through enhanced am- 
bipolar diffusion, which produces filamentary supercritical structures that 
are generally long lived. Only the densest parts of the filaments — the dense 
cores — are directly involved in star formation. Even the dense cores are sig- 
nificantly magnetized, with flux-to-mass ratios about half the critical values 
or more. The strong magnetization ensures that ambipolar diffusion con- 
tinues to play a role in the core evolution, reducing t heir iri t ernal velocity 
dispersions to typically subsonic levels, as observed ( Mvers . 19951 ). Since 
only a fraction of the cloud material is magnetically supercritical (and thus 
capable of forming stars in the first place) , and only a small fraction of the 
supercritical material is directly involved in star formation, the efficiency of 
star formation is naturally low. 

An attractive feature of the "turbulence-accelerated, magnetically regu- 
lated" star formation in diffuse, subcritical clouds (or sub-regions) is that 
the stars are expected to form at distributed locations. This is because 
the self-gravity, although important on the scale of individual supercritical 
cores, is canceled out to a large extent, if not completely, by magnetic forces 
on large scales. The cancelation makes the clustering of dense cores difficult, 
unless the cores happen to be created close together by converging flows. By 
the same token, star clusters are more likely formed in dense, self-gravitating 
clumps that are magnetically supercritical as a whole, as we discuss below. 

Even though the principle of the "turbulence-accelerated, magnetically 
regulated" star formation is straightforward, and well illustrated by 2D sim- 
ulations, much work remains to be done to flrm it up. Among the needed 
refinements is the extension o f the 2D calcu l ation s to 3D. A first step in this 
direction has been taken by iKudoh et al.l (|2007l ) and Nakamura & Li (in 
preparation). Another improvement would be in the maintenance of turbu- 
lence, which will be more difficult, given our limited understanding of the 
origins of turbulence in molecular clouds in general, and in diffuse regions 
of distributed star formation in particular. 

A region of inefficient, dis tributed sta r form ation of considerable current 
interest is the Pipe nebula. Alves et al. ( 200?! ) found more than 150 dense 
quiescent cores with a mass distribution that resembles the stellar IMF. The 
cores are distributed along a filament of more than 10 pc in length. None of 
the cores except the most massive one has yet collapsed and formed stars. 
This is difficult to understand, given their rather short free fall times, unless 
the core creation is well synchronized. The synchronization is most naturally 
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done by a large scale shock. In such a case, the strong compression is ex- 
pected to induce fast motions inside the filaments and cores, unless the cores 
are magnetically cushioned. The Pipe nebula may turn out to be a good 
example of compression-induced, magnetically regulated star formation. 



1.3.2 Turbulence, Gravity, and Cluster Formation 

The most detailed simulations of cluster formation are performed using the 
SPH technique, which gene rally does not include magnetic fields (see, how- 



ever. 



Price and Batd . 120071 ) . These simulations are well suited for studying 



the interaction between turbulence compression, gravi tational collapse, and 



especially mass accretion onto collapsed objects (see iBonnell et al.l . 12007 



review). They fall into two categc 


Bate et al. 




2003|; 


Bonnell et al.. 



turbulence ([Klessenl . l200ll ). In the former, the initial turbulen ce typically 



gener ates several den se clumps, in which small stellar groups (jBate et al 
2003h or sub-clusters dSonnell et al.l . bood ) form, depending on the number 



of Jeans masses initially contained in the cloud. The most massive member 
of each sub-cluster typically sits near the bottom of the gravitational poten- 
tial well of the local collection of gas and stars, gaining mass preferentially 
through competitive accretion. The sub-clusters of stars and gas merger to- 
gether at later times, making deeper gravitational potentials where the more 
massive stars near the bottom of the potential well can accrete preferentially, 
growing to higher masses. 

The most attractive feature of the non-magnetic SPH calculations of clus- 
ter formation is that an IMF that broadly resembles the observed one is 
produced. In particular, the characteristic mass near the knee of the IMF 
is attributed to the initial cloud Jeans mass (jBonnell et al.l . hood ). and the 



steeper, Salpeter-like mass distribution above the knee is mostly shaped by 
competitive accretion. A potential difficulty is that the Jeans mass (and 
thus the characteristic mass of IMF) is expected to vary from region to re- 
gion because of variations in both density and temperature. There is little 
observational evidence to support this expectation. The difficulty can be 
alleviated to some extent by the heating and cooling behaviors of molec- 
ular gas, which tend to set a characteristic density below (above) which 
the temperature de creases (increases) with density (jJappsen et al.l . l2005l : 



Bonnell et aP - bood ). Whether the thermal regulation of Jeans mass works 



for very dense clusters wh ere the average Jeans mass is expected to be very 
small remains to be seen (jMcKee and Ostrikerl . [20071 ) . 

A problem with cluster formation in decaying turbulence is that star for- 
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mation may be too rapid. According to ( Klessen . 200ll ) ~ 20 — 30% of the 
cloud mass is accreted by sink particles in one free fall time tfj. The value of 
tff depends on the density (and thus the initial Jeans mass Mj for 10 K gas). 
For Mj = 1 Mq (needed to reproduce the knee of IMF), t// ~ 10^ yrs. This 
may be short compared to the lifetimes (typically 1 Myrs or more) estimated 
for nearby embedded clusters, particularly those that include a substantial 
population of relatively evolved Class II and Class III objects. Increasing 
the free-fall time to 1 Myrs would increase the Jeans mass to 10 Mq, which 
would produce too high a mass for the knee of the IMF, unless the conversion 



of dense, self-gravitating gas int o individual stars is inefficient (jShu et al 



19991 : iMatzner and McKeel . bood ^ . Alternatively, the rate of star formation 



can be reduce d significantly, if the turbulence is constantly driven on a small 
enough scale ( Klessenl . l200ll ). The small-scale driving may, however, modify 
the hierarchical cluster formation observed in the simulations of the decaying 
turbulence case, an d affect compet itive accretion and thus the high mas s end 
of the IMF. Indeed, Klessen ( 200ll ) and Vazquez-Semadeni et al. ( 20031 ) sug- 
gested that turbulence driven on small scales may correspond to distributed 
star formation, whereas that on large scales to clustered star formation. 
According to these works, large-scale driving does not slow down star for- 
mation significantly, so the problem of too rapid star formation in clusters 
may remain (but see Sect. 11.41 below). The problem ca n be alleviated some- 



what by the inclusion of a strong magnetic field (e.g.. iHeitsch et al.l . 12001 



Vazquez-Semadeni et al. . 20051 : Tillev and Pudritz . 2007 ). In the next sub 



section, we examine the role of protostellar outflows in slowing star for- 
mation by replenishing the turbulence in localized regions of active cluster 
formation. 



1.3.3 Cluster Formation in Protostellar Turbulence 

The possibility of outflows replenishing the energy and momentum dissi 



pated in a star-forming cloud was first examined in detail bv lNorman and Silk 



( 1980l ). They envisioned the star-forming clouds to be constantly stirred up 
by the winds of optically revealed T Tauri stars. The id[ ea was strength- 
ened by the discovery of molecular outflows ( McKee . 19891 ). which point to 
even more powerful outflows from the stellar vicin i ty during the embedded , 



protostellar pha se of star formation ( Lada . 1985 : Bontemps et al. . 19961 ) 



Shu et al.l (|l999l l estimated the momentum output from protostellar out- 
flows based on the Galactic star formation rate, and concluded that it is 
sufficient to sustain a level of turbulence of ~ 1 — 2 km/s, similar to the 
line widths observed in typical GMCs. If the majority of stars are formed 
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in localized pa rsec-scale dense cl umps that occupy a small fraction of the 



GMC volume (jLada et al.l . Il99ll ). their ability to influence the dynamics 
of the bulk of the GMC material will probably be reduced; other means of 
turbulence maintenance may be needed in regions of relatively little star for- 



mation, as concluded bv I Walawender et alj (|2005l ) in the case of the Perseus 



molecular cloud. The concentration of star formation should, however, make 
the outflows more important in the spatially limited, but arguably the most 
interesting regions of a GMC — the regions of cluster formation, where the 
majority of stars are thought to form. 

The importance of outflows on cluster formation can be illustrated using 
a simple estimate. Let the masses of the cluster-forming dense clump and 
the stars formed in it be Mc and M^,, respectively. If the outflow momentum 
per unit stellar mass is P*, then there is enough momentum to move all of 
the clump mass to a speed 

,P. = ,kn./s(^](,^]. (1.1) 



Mc ' VO.l/ V50 km/s 

where e = M^/Mc is the star formatio n efficiency (SFE). Fo r embedded 



clusters, the SFE can be 10% or more (jLada and Ladal. 120031'). Th e value 



of P=K is somewhat uncertain. For low mass stars, iNakamura and Lil (j2007l ) 



estimated a plausible range between 10 — 100 km/s. For the best studied 
CO outflow-driving early B stars, the value is estimated at 50 — 100 km/s 
(D. Shepherd, priv. comm.). For the fiducial values of e and P* adopted in 
equation (jl.ip . the outflow-driven average velocity is comfortably above the 
typical turbulence velocity of 1-2 km/s observed in nearby cluster-forming 
clumps. Indeed, if all of the cluster members were to form simultaneously, 
the clump would quickly become unbound. If the stars are formed more 
gradually (as evidenced by the presence of objects in a wide range of evolu- 
tionary states, from prestellar cores to Class III sources), then there is the 
possibility for the outflows to replenish the dissipated turbulence, keeping 
the cluster formation going for a time longer than the global free fall time. 

Quasi-equilibrium cluster formation in outflow-driven turbulence is illus- 
trated in Fig. 11.61 It is clear that dense cores tend to collect near the bottom 
of the gravitational potential well, where most of the stars form. Most of 
the outflow momentum is coupled, however, into the envelope, where most 
of the clump mass resides. The coupling of momentum into the envelope 
is facilitated by outflow collimation, which enables the outflows to propa- 
gate further away from the central region, and drive motions on a larger 
scale that would decay more slowly. Gravity also plays an important role 
in generating the turbulent motions by pulling the slowed down material 
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Fig. 1.6. Snapshots of a 3D MHD simulations of cluster formation including feed- 
back from protostellar outflows. Plotted are the volume density map, velocity 
vectors, and contou rs of gravitational potent ial on a slice through the simulation 
box. For details, see iNakamura and Lil (j2007t ). 



towards the center and setting up a mass circulation between the core and 
envelope. For parameters appropriate for nearby embedded clusters such 
as NGC 1333 (~ 10^ Mq in mass and ~ 1 pc in size), the gravitationally 
induced infall and outflow-driven expansion are roughly balanced, creating a 
quasi-static environment in which the clump material is converted into stars 
at a relatively leisure pace. For the particular example shown, the average 
star formation efficiency per global fre e fall time is abou t 3%. This value 
is in the range inferred for NGC 1333 (jNakamura and Lil . HqO^) and other 
objects (jKrumholz and Tanl . 120071 ) . 

The protostellar turbulence in regions of active cluster formation is a spe- 
cial type of driven turbulence. It is driven anisotropically over a range of 
scales, as the collii nated outflow propagates progressively further away from 
the driving source (jCunningham et al.l . l200^ : lshang et all . bood : | iBanerjee et al 
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20071 ). There is a break in the energy spectrum n ear the characte ristic length 



scale of the outflows, as predicted analytically bv lMatznerl (j2007l ). The break 
may provide a way to distinguish the protostellar turbulence from other 
types of turbulence. 

It is interesting to speculate how massive stars may fo rm in protostellar 



outflo w-driven turbulence. One possibility, pointed out in lLi and Nakamura 



(|2006l l. is that, when the central region of the cluster-forming clump becomes 
dense enough, the outflows from the low mass stars may become trapped. 
The trapping of outflows reduces the stellar feedback into the massive enve- 
lope, which leads to more mass falling to the central region, which in turn 
makes it more difficult for the outflow to get out. It may lead to a run-away 
collapse of the central reg ion, with a large i nass infall rate that is conducive 
to massive star formation (|McKee and Tanl . |2002| ) . The massive stars formed 
in this scenario will naturally be located near the bottom of the gravi tational 
potential well, for wh ich there is now growing observational evidence (jGaravl . 
2005l : ICesaronil . l2005l ^. The trapping of outflows as a cluster-forming clump 
condenses may also provide a natural explanat ion for the accelerating star 
formati on inferred for the Orion Nebula cluster ( Huff and Stahled . 120061 ) and 
others (jPalla and Stahleil . bood ^ . Detailed numerical simulations are needed 
to firm up this supposition. 



1.4 Global Theory of Star Formation in Turbulent Clouds 

Giant molecular clouds and their smaller scale constituents, star forming 
clumps and prestellar cores, obey several remarkable scaling relations, which 
indicate that star formation to a significant extent exhibits a generic behav- 
ior. The most important of these scaling relations are: 



(i) 



(ii) 
(iii) 



Larson's velocity-size relation ( Larsonl . Il979l ). which expresses that 
the (3D) rms velocity disper sion cr^.^n of s tructures o f size L scales as 
Oj2^^jn^ L^", where ^ 0.4 ( Larsonl . Il98ll ). Recently Hever and Brunt 



(|2007l ) have shown that the relation is very tight when measured in 
a consistent manner over a limited range of scales, and that the nor- 
malization is essentially the same in violently s tar-forming and nearly 



quiescent molecular clouds ( Hever et al. . 20061 ) . 



Larson's velocity- mass relations (jLarsonl . Il98ll ). which expresses that 
the densest structures on each scale are in near-virial equilibrium. 
The cluster mass functio n, which is a power la w with slope of dN/dM) 
slightly less steep that -2 (jDowell et al.l . l2006l . and references therein) — 
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so stellar birth is distributed nearly evenly over mass, although a 
slightly larger number of stars are born in rich clusters, 
(iv) The Core Mass Function (CMF) and the Initial Mass Function (IMF), 
which are approximate power laws for large enough masses. Individ- 
ual examples may show deviations, due to statistics or other influ- 
ences, but the general power law trends are well established. 

These structures thus appear to be members of a hierarchy of scales, char- 
acterized by approximate scaling laws. Two additional, and crucially im- 
portant, properties of this hierarchy are that: 

(a) On each scale the actual lifetimes of the structures appear to be 
longer than their estimated turbulent decay times (as given by the 
turbulence c rossin g times). To maintain their (rather precise; cf. 
Hever et all , bood l velocity normalization in the hierarchy, the mo- 



tions on each scale thus need significant driving, 
(b) A related but largely independent property is that the mass drain 
due to the star formation process is approximately constant across 
scales, at a level o f only a few percent per unit free fall time (cf. 



Krumholz and Tanl . 120071 . and references therein) 



The various properties enumerated above are linked in ways that are likely 
not to be coincidences: items (a) and (b) above imply that the draining of 
mass from each scale per unit crossing time is quite small. This property 
is in principle independent of the property that the structures appear to be 
long lived; one could imagine the structures to have a small draining rate but 
still be quite short-lived. The opposite could also be the case; the structures 
could be drained more rapidly by star formation but could still be long lived 
if the geometrical shapes of structures were somehow maintained as they 
were being drained. 



1.4-1 Questions and Difficulties 

Most of these propertie s are difficult to exp l ain: t he small rate of star forma- 
tion per free fall time ( Krumholz and Tan . 20071 ). the driving source of the 
motions, the regularities of the cluster mass function, the core mass function, 
and the stellar initial mass function are all long standing problems. 

One can attempt to "explain" Larson's density and velocity relations as a 
consequence of "virialization" , since they are roughly consistent with virial 
equilibrium. With velocity scaling as L^^/^ (rather than the observed ex- 
ponent ~ 0.4) and density as L~^, one could imagine that gravity main- 
tains the velocity spectrum, as structures shrink in size and descend into 
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ever deeper potential wells. That, however, still begs the question "why 
would that happen" ; specifically, why would this happen in such a way that 
the life times remain long and the draining rates at each level in the hi- 
erarchy remain nearly constant (and low)? One would rather expect that 
the characteristic time scale of such a cascade would be similar to the local 
gravitational free fall time. The driving- by-gravity explanation is also very 
difficult to reconcile with the fact that potential wells actually get shallower 
with decreasing size, rather than deeper! 

Likewise, it is very difficult to explain how local feedback could, without 
exceptional fine tuning, both maintain a power law scaling of the velocity 
with size — the energy input rate would have to match the turbulent dis- 
sipation rate and its variation with scale over a large range of scales — and 
maintain the proper amount of support against gravitational collapse to keep 
the draining rate constant and low over a range of structure sizes. 

One is lead to conjecture then that local feedback, rather than maintaining 
the observed scaling laws over a range of scales instead may be responsible 
for breaking those scaling laws at lo cal sites with rapid s tar fo rmation; as 
observed in for example CS clumps ( Krumholz and Tan . 200?! ). where the 
mass density (and possibly also draining rates) exceed values expected from 
the scaling relations. But that still leaves the questions of explaining the 
power law scalings that hold at most scales and locations, of explaining 
the maintenance of these motions against dissipation, and of explaining the 
constancy of draining rates over all other scales. 

In fact, most of the scaling and lifetime properties may be explained, in 
a consistent and complete manner, by essentially "turning the explanation 
up-side-down"; it is not the gravity that is in the lead but the kinetic en- 
ergy, cascading from large to small scales, with gravity picking up those 
exceptional regions of space where the local gravitational potential energy 
dominates the local kinetic energy (i.e. the local virial parameter, as given 
in Eq. 11.21 happens to be small — cf. Fig. II. 7p . 



1.4 -2 A Turbulent Cascade Origin of the Scaling 

A suitable entry point for discussion of the variou s scali ng p roperties is Lar- 



son's v elocity-size relation. As shown by Boldvrev ( 20021 ) and Boldvrev et al 



(I2OO2I I such a scaling follows inevitably from the scaling relations of super 
sonic turbulence. A number of observationally well defined diagnostics are in 
agreement with the corresponding diagnosti cs extracted from simulations of 



super sonic (and super- Alfvenic) turbulence (IPadoan et al .1. 1 1 99 71. 1 1 9981. |200 1 



20041 ). One can even measure, using the method of iLazarian and Pogosvan 
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the velocity power spectrum exponent of the turbulence; the result is 
a wavenumber expone nt ~ —1.8, coiisistent with a velocity-size scaling with 
an exponent ^ 0.4 (jPadoan et al.1 . bood l . 



So, what we are observing as Larson's velocity-size scaling law is, most 
likely, just the inevitable and robust consequence of a cascade of kinetic 
energy from the largest (injection) scales, which occur at scales of the order 
of the galactic disk thickness, where kinetic energy may be contribu t ed by 
several sources (de Avillez and Breitschwerdt . 2005. 2007: Kim et al.. 2006: 



Qstriked . l2007l : iMcKee and Ostrikeil . l2007l l. No matter how that energy 
is fed into the medium it must cascade to very small scales before it is 
dissipated. This occurs, in supersonic turbulence as well as in subsonic 
turbulence, via a turbulent cascade across an "inertial range", where the 
denomination "inertial" signifies that the motions are maintained by inertia. 

This explains both the power law distribution of velocity (and hence Lar- 
son's velocity-size relation), and th e source of ap parent "driving" over a 
range of scales. As already shown bv iLarsonl ( 19791 . cf. his Fig. 1) the power 
law scaling continues to sizes comparable to the thickness of the galactic disk, 
where it exhibits a break (at larger scales the velocity dispersion is mainly 
associated with fluctuations in the galactic disk rotation rate). The title of 
Larson's (1981) paper was indeed, and very appropriately, "Turbulence and 
Star Formation in Molecular Clouds"! 

The super-Alfvenic turbulence sce nario may also be used to de rive the 
form of the stellar IMF analytically ( Padoan and Nordlund . 20021 ) and to 
explain how brown dwarfs can form by the same mechanism as normal stars 
(jPadoan and Nordlund . |2004| ). 



1.4-3 Interaction of Turbulence and Self gravity 

One might object to the comparisons of diagnostics from simulations of 
supersonic turbulence without selfgravity with observations of molecular 
clouds since these, according to the previous discussion, are often close to 
virial equilibrium, and it would thus seem that selfgravity cannot be ignored. 

One could, on the other hand, also argue that the success of such ex- 
ercises might indicate that selfgravity is only important in local regions, 
and that these are sufficiently small not to disturb the comparisons signif- 
icantly. Such a conjecture can actually be confirmed from an analysis of 
loc ally evaluated virial rat i os ('v irial parameters'), which we may define as 
(cf. iKrumholz and McKeel . boosl ) 

^'^rms /I o^ 
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Fig. 1.7. The mass weighted probabihty distribution of a (left), normalized to the 
virial ratio for the whole domain, in a numerical simulation of driven Mach=18 
supersonic turbulence with selfgravity (Padoan & Nordlund in prep.; resolution 
1000^), and the dependence of the star formation rate on the large scale a (right) for 
a number of such simulations (stars — the broken power law curve is only suggestive). 

where Vxms is a local (3-D) rms velocity dispersion, measured in some neigh- 
borhood of size L around comoving, lagrangian tracer points where the local 
average density is p. 

The left hand side panel in Fig. 1 1.71 shows an example of the mass weighted 
probability distribution of the local virial ratio, using values of v'^^ns ^^'^ P 
evaluated locally. It is obvious that a large fraction of the mass has local 
values of a much larger than the virial ratio ao based on the rms velocity 
and average density for the whole domain, and that only a small fraction of 
the mass (~1% in this case) has a smaller value. 

A consideration of how local shearing motions are modified by compression 
explains this behavior. Suppose we have identified a region of size L with a 
small density enhancement (and thus potentially a selfgravitating 'object'), 
which initially has a value of a close to the value for the whole box and a 
density close to the mean density. If we choose L to be for example 1/4 of 
the size of the full box then v must likewise be 1/4 of the average rms value 
to obtain the same a, so the initial selection is already limited to a subset 
of the volume (in the rest of the volume a local a evaluated in this fashion 
would be larger, even initially). 

In a local neighborhood, and for a limited time, one can assume the local, 
comoving dynamics to be subsonic. As can be confirmed by data-browsing, 
the velocity field in supersonic turbulence consists of patches with smooth 
flow, separated (ideally) by discontinuities of the velocity and the mass den- 
sity. With smooth velocity across a patch one can subtract off the mean 
speed and consider the local velocity field, which can be split into a solenio- 
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dal and a compressive part. In a sufficiently small neighborhood in space 
the flow remains subsonic with respect to the mean, and in a sufficiently 
small neighborhood in time — until the patch hits a shock — one may thus 
consider the local dynamics to be subsonic, smooth and continuous. 

One consequence of this is that, within a thus limited neighborhood in 
space and time it makes sense to think in terms of approximate conservation 
of local angular momentum. The subsonic nature of the local dynamics 
means, for example, that convergent motions couple to shear and rotations 
in the well known sense that contraction tends to cause spin-up. 

Now suppose the region in question is compressed to a tenth of its lin- 
ear size, with compression taking place in two spatial directions. The cross 
section decreases with a factor of 100, but the density increases correspond- 
ingly, so the denominator in Eq. I1.2l is unchanged. Because of the tendency 
to locally conserve angular momentum the v'^ in the numerator tends to 
increase by a large factor, and one flnds that the a characterizing this local 
region has grown tremendously, which prevents the object from becoming 
selfgravitating. 

In general one finds that in regions of compressions, which are otherwise 
the locations that are most favorable for creating large mass densities, the 
tendency for growth of vorticity more than counters the growing influence of 
selfgravity, which makes it hard to achieve collapse under selfgravity, even 
when the virial parameter measured on the large scale is of the order unity. 

There is thus a purely dynamical mechanism available to explain why 
structures on all scales in GMCs tend to survive longer than would be in- 
dicated by estimates based on the large scale virial parameter (such an 
estimate is exactly what is being used when one says "it is remarkable that 
these structures live longer than a few free fall times"). The mechanism is 
even unavoidable; there is no way that this could not happen, in cases where 
the virial parameter based on the large scale is larger than unity! 

When, on the other hand, the large scale virial parameter is less than unity 
opportunities rapidly open up for local collapse; the fraction of compressive, 
converging flows that already initially have low a increases, and a rapidly 
growing number of compressed regions flnd themselves to have local virial 
ratios low enough for collapse. To a flrst approximation one can see the 
dependence of the draining rate on a (right hand side panel of Fig. II. 7p 
as resulting from a shift of the PDF of a (left hand side panel of Fig. 11.71 
considered as a function a rather than a/ao). 

Such an argument may be applied recursively. Structures that are much 
smaller than the initial size but that still have low enough local virial ratios 
for collapse can be taken as a basis for 'renormalization'. One concludes that, 
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for any large scale virial ratio, a hierarchy of smaller structures exists, with 
local virial ratios of the order unity, and with internal structures consisting 
of yet smaller structures with local virial ratios of the order of unity. The 
hierarchy must be such as to conserve the draining rate at each scale; since 
the smaller scale free fall times are shorter than the larger scale free fall 
times the smaller scales are "slaved" to the larger scales. 



1.4 -4 Self regulated Star Formation 

One very important aspect of the difficulty to achieve locally small virial 
ratios when the large scale virial ratio is larger than unity is that the rate 
of gravitational collapse of mass out of this hierarchy (which presumably is 
a good proxy of the star formation rate — here denoted SFR), turns out to 
be extremely sensitive to the large scale virial ratio. 

The right hand side panel in Fig. 11.71 illustrates the point, by showing a 
summary of the star formation rates achieved in a number of simulations 
where the large scale alpha was varied, both by changing the mean density 
and by changing the Mach number (Padoan & Nordlund, in prep.). The 
functional behavior is such that the SFR is approximately proportional to 
for a larger than unity, while for smaller alpha the SFR is roughly 
c onstant and very large 30 % per free fall time). 

Krumholz and McKeel (I2OO5I I arrived at an analytical estimate of the de- 



pendence of the SFR on a that is somewhat less ste ep but qualitat ively 
similar, based on results f rom numer ical simulations by iKlessen (l200lh and 
Vazauez-Semadeni et al.1 (|2003l . I2OO5I ). 



In the regime where the SFR is nearly constant it is si gnificantly larger 
than values characterizing star formation in the Milky Way ([Krumholz and Tanl . 



20071 ). We thus conclude that galactic star formation is operating in a regime 



where a is on the steep, power-law part of the dependence. 

Does that then not require an unlikely fine tuning? Why would the large 
scale a happen to have a value such that the SFR agrees with the observed 
one? Well, suppose that it did not: if the value of a were too large there 
would be practically no star formation, the driving of the ISM from massive 
stars and SNe would cease, the level of turbulence would be reduced, the 
scale height of the ISM would drop and the large scale a would decrease — 
both because of reduced velocities and because of increased densities. 

Conversely, if a were too small the SFR would be too high (relative to 
the level required to sustain the ISM turbulence), and the situation would 
rectify itself via the same feedback mechanism. 

We conclude that the mechanism that determines the level of the star for- 



1.5 Concluding Remarks 



23 



mation rate is large scale stellar feedback. Note also that this only works if 
the driving of the interstellar medium turbulence is, to a significant but not 
necessarily exclusive extent, due to large scale feedback from star formation. 
That this is so is already abundantly established by elaborate numerica l 
simu lations, in particular the ones by Ide Avillez and Breitschwerdtl (|2005l . 
20071 ) . Such a scenario is also entirely consistent with the threshold assump- 



Conrov et al. 


2007^; 


Naeai et al. 


( Kennicuttl, 


1998|, 


2007|). 



1.5 Concluding Remarks 

The last decade has brought significant observational and theoretical progress 
to the study of star formation in molecular clouds, particularly in observa- 
tions of embedded clusters and simulations of turbulent clouds. With large 
scale surveys of both molecular clouds and their populations of protostars 
and pre-main sequence stars now in hand, and with increasingly sophisti- 
cated numerical models, there is excellent potential to rapidly advance our 
understanding of the macrophysics of star formation. 

To do so will require better means for comparing observations and simula- 
tions. Observers and modelers need to devise a set of common "observable" 
statistics that may be used to perform quantitative comparisons. Given the 
advances in measuring the distribution, evolutionary state, and even kine- 
matics of YSOs and dense cores, it is of great interest to relate the properties 
of the sink particles in simulations of star formation to the observed prop- 
erties of dense cores and YSOs. In particular, it is important to determine 
whether the observed spatial distributions of protostars and more evolved 
pre-main sequence stars can be reproduced. The evolution of clusters should 
also be explored through N-body simulations with realistic and time variable 
gas potentials to better understand the dynamical states implied by their 
observed morphologies. 

There are areas where significant progress is needed on the observational 
front. Much progress is still needed to understand the ages and age spreads 
of embedded populations of YSO. These are an important constraint on 
the lifetimes of clouds and the duration of star formation. Significant ad- 
vances have been made in the measurements of magnetic fields; however, 
since these measurements often target regions of active star formation, the 
field structure in the bulk of cloud material is poorly constrained. Large- 
scale mapping of the strength and direction of the magnetic field is needed 
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to ascertain its dynamical role, particularly in relatively diffuse regions of 
clouds, where the field may play a dominant role. 

Theoretically, there is a need to tie local simulations of star formation to 
global simulations of structure formation on GMC scales and beyond. On 
the GMC and smaller scales there are additional effects not included in this 
review. These include the generation of sma, ll scale structure from global 
cloud collapse ( Burkert and Hartmann . 20041 ) . generation of structure by 



thermal instabilities at the interfaces between the warm ISM and molecular 
clouds (Heitsch et al., 2007), self-shielding and the transition from atomic 
to molecular phase in the interior of molecular clouds. In order to avoid 
the need for ad hoc initial and boundary conditions such modeling may be 
performed as sub-sets of larger scale and longer duration global simulations, 
as often done in the study of galaxy formation in the context of cosmological 
structure formation simulations. 

On larger scales, simulations need to account for the life cycle of molec- 
ular gas, and determine whether "lifetime" is even a meaningful concept in 
the context of molecular clouds. Relevant questions to investigate include 
"Given a snapshot of a molecular cloud, how long time does it take before 
half of the molecular gas present in the snapshot is recycled?" , "How much 
of that went into stars?" , "How much was turned into warm and hot ISM, 
respectively" , "How much new molecular gas was added in the mean time?" , 
"How much of the feedback is due to winds, UV- radiation, and supernovae 
and on what scales is this feedback deposited?" . Such questions may be an- 
sw ered by realistic high resolution glob al ISM simulations (such as the ones 
by de Avillez and Breitschwerdt . 20071 ). which include the vertical structure 
of the galactic disk, realistic cooling functions and UV heating, selfgravity, 
and magnetic fields. A next step in such simulations could be to let the star 
formation rate be determined self-consistently, by including approximations 
of how the star formation rate depends on the virial ratios of star forming 
structures and their average magnetic fields. 
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